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ABSTRACT 


The  CTBT  negotiations  have  endowed  small  seismo/acoustic  events  (mb  ~  2.5)  with  a  much  greater 
significance  than  they  have  had  in  the  past  and  have  increased  the  need  for  automated  regional  dis¬ 
crimination.  We  have  been  contracted  by  AFTAC  to  develop  a  small  event  discriminant  that  uses  a 
time-frequency  expansion  (sonogram)  of  seismic  coda  to  discriminate  ripple-fired  mining  explosions 
from  single  explosions  and  earthquakes.  The  Automated  Time-Frequency  Discriminant  (ATFD)  which 
we  have  developed  uses  a  binary  sonogram  which  is  derived  from  the  original,  spectral,  sonogram  by 
the  application  of  filters  which  replace  the  spectral  information  with  a  binary  code  which  simply  re¬ 
flects  local  spectral  highs  and  lows.  We  have  found  that  the  binary  patterns  which  we  extract  from  the 
coda  of  ripple-fired  events  are  banded  and  thus  distinct  from  those  obtained  from  single  explosions  and 
earthquakes.  The  bands,  which  result  from  source  finiteness,  intershot  delays  or  a  combination  of  the 
two,  are  largely  independent  of  time  and  of  the  recording  component.  The  ATFD  uses  three  statistical 
tests  to  measure  the  time-  and  recording  component-independence  and  automatically  recognize  these 
bands.  All  of  the  raw  discrimination  parameters  produced  by  these  tests  are  merged  into  a  single 
discriminant  score  with  multivariate  statistics. 

During  the  past  year  we  have  tested  the  ATFD  on  data  collected  by  single  stations,  arrays  and  regional 
networks  in  central  Asia,  Europe  and  North  America.  We  have  found  the  automated  discriminant  to  be 
robust,  with  misclassification  probabilities  ranging  from  0.5  to  3.5  percent.  The  method  is  easily 
adapted  to  any  kind  of  regional  seismic  data.  No  expert  analysis  of  the  time-frequency  patterns  is 
required. 

To  improve  the  ATFD  we  have  analyzed  wavelet  expansions  of  coda.  We  have  found  that  wavelets  are 
capable  of  much  higher  resolution  in  time  than  is  possible  with  a  multitaper-sonogram  expansion.  We 
are  currently  using  a  new  bias-minimizing  wavelet  technique  (Lilly  &  Park,  1995)  to  process  3-compo¬ 
nent  data  and  analyze  the  evolution  of  spectral  amplitude  and  polarization  with  time  and  frequency  and 
improve  our  understanding  of  ripple-fired  coda.  We  intend  to  determine  if  a  low  frequency  spectral 
signature  (perhaps  caused  by  source  fmiteness)  might  be  used  for  discrimination  at  far-regional  dis¬ 
tances.  We  have  conducted  a  preliminary  survey  of  discrimination  at  low  frequencies  using  the  Kazakh 
NRDC  network. 

We  have  begun  to  look  at  discriminants  based  on  hybrid  (seismo/acoustic)  data  and  will  shortly  collect 
data  to  help  us  address  the  issue  of  outliers.  These  data  will  be  collected  this  summer  in  Wyoming  and 
will  consist  of  broadband  seismic,  acoustic  and  infrasound  recordings  of  ground  truthed  mining  explo¬ 
sions  in  the  Black  Thunder  coal  mine.  With  this  experiment  we  will  also  be  examining  the  utility  of 
continuous  GPS  for  event  detection.  At  regional  distances  from  the  Black  Thunder  mine  we  will 
deploy  a  network  of  3  continuous  GPS  receivers  and  will  use  the  data  to  look  for  ionospheric  perturba¬ 
tions  resulting  from  the  near- vertical  acoustic  emissions. 

Keywords:  comprehensive  test  ban  treaty  (CTBT),  time-frequency  expansion,  wavelets,  automated 
discrimination,  multivariate  statistics,  coda  cepstrum,  ripple-firing,  binary  sonogram. 
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A.  OB.TECT1VES/PROGRESS 


(1)  Objectives 

During  the  term  of  our  existing  contract  we  will  expand  significantly  research  we  have  conducted  into 
the  small-event  discrimination  problem.  Our  interest  has  been  in  using  the  spectral  characteristics  of 
whole  regional  waveforms  to  discriminate  ripple-fired  from  non-ripple-fired  events  (inch  earthquakes 
and  single-event  explosions).  Our  objectives  are  to  apply  our  automated  whole  waveform  time-fre¬ 
quency  discriminant  (ATFD)  to  large  event  populations  in  varied  data  sets  to  test  transportability  and 
robustness  and  to  enhance  the  ATFD  with  more  sophisticated  processing.  In  specific  we  have  the 
following  objectives: 

LI  Robustness  and  transportability 

Apply  our  existing  ATFD  (Hedlin  et  al.,  1990)  to  a  number  of  dissimilar,  well  separated,  re¬ 
gional  data  sets  with  large  populations  of  ripple  and  non-ripple-fired  events  to  gauge  robustness 
and  regional  dependence.  We  will  apply  the  technique  to  vertical  component  data  from  several 
networks  and  arrays  using  low  (1  to  20  Hz)  and  broad  (1  to  100  Hz)  frequency  bands.  Particu¬ 
lar  attention  will  be  paid  to  the  cause  of  outliers.  We  will  determine  the  extent  to  which  array 
data  can  suppress  noise  and  increase  the  range  of  the  ATFD. 

1.2  Software  development 

Enhance  the  ATFD  via  (1)  wavelet  analysis,  (2)  more  advanced  spectral  analysis  techniques 
[e.g.  the  statistics  of  Higher  Order  Crossings  (HOC)],  (3)  advanced  processing  techniques  to 
permit  full  use  of  modem  three-component  (3C)  networks  and  arrays.  3C  data  sets  examined 
under  objective  1.1  will  be  re-analyzed  to  assess  improvement. 

1.3  Comparisons  with  other  techniques 

Analyze  the  same  data  sets,  discussed  under  objective  1.1,  with  a  complementary  technique 
(e.g.  a  “regionally  trained”  spectral  ratio  method)  to  assess  relative  capabilities  under  different 
settings.  Our  intent  is  to  develop  a  complementary  technique  that  might  be  merged  through 
evolutionary  programming  with  our  own  to  provide  a  more  comprehensive,  multivariate,  dis¬ 
criminant.  The  time-frequency  approach,  taken  alone,  will  not  discriminate  between  single¬ 
event  explosions  and  earthquakes. 

1.4  Low-frequency  discrimination 

Address  the  question  of  what  low  frequency  (e.g.  1  to  20  Hz)  time-independent  signature 
should  be  produced  by  mines  that  use  short  (e.g.  20  ms)  delays  especially  when  delay  times  are 
irregular.  Use  any  available  “ground  trath”  data  to  explain  any  modulations  observed  at  low 
frequencies  and  test  the  theory  that  they  might  be  due  to  temporal  finiteness  of  the  source 
(Hedlin  et  al.,  1990). 

1.5  Contribution  of  software  to  database  accessing  systems 

Develop  algorithms  that  are  designed  to  operate  on  the  CSS  3.0  database  structure  and  will  be 
available  to  all  interested  parties. 
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(2)  Progress 


The  automated  time-frequency  discriminant.  Under  previous  Air  Force  contracts  (F19628-87-K-0013 
and  F19628-88-K-0044)  we  developed  a  discriminant  that  seeks  long-lived  spectral  modulations  in 
major  seismic  phases  and  coda.  Long  lived  modulations  can  be  produced  by  seismic  resonance  and  by 
ripple  firing.  We  expand  time  series  into  time-frequency  displays  (sonograms)  to  search  for  these 
patterns.  Hedlin  et  al  ( 1989)  developed  a  procedure  whereby  a  binary  sonogram  is  derived  from  the 
original,  spectral,  sonogram  by  the  application  of  filters  which  replace  spectral  information  with  a 
binary  code  which  simply  reflects  local  spectral  highs  and  lows.  Hedlin  et  al.  (1990)  produced  a 
procedure  to  automatically  recognize  time  independent  patterns.  This  technique  utilizes  a  two-dimen¬ 
sional  Fourier  transform  of  the  binary  sonogram  which  reveals  the  dependence  of  the  binary  pattern  on 
frequency  and  time.  In  view  of  its  resemblance  to  the  cepstrum  (which  identifies  periodicities  in  single 
spectra),  and  the  fact  that  it  is  derived  from  onset  and  coda  phases  we  now  refer  to  it  as  the  coda 
cepstrum  (Hedlin  et  al.,  1995). 

Hedlin  et  al.  (1989)  also  noted  that  the  modulations  were  remarkably  independent  of  the  recording 
direction.  As  an  example  of  this  we  display,  in  figures  1, 2,  and  3,  a  summary  of  an  analysis  of  record¬ 
ings  of  earthquakes  and  quarry  blasts  made  by  the  KNET  in  Kyrgyzstan.  Examining  the  quarry  blasts 
(located  just  northeast  of  CHM)  and  the  earthquakes  located  to  the  southwest  of  the  network  we  found 
that  the  network  average  cross  correlations  of  the  binary  patterns  between  the  recording  components 
(average  of  EW  vs  Z,  Z  vs  NS  and  EW  vs  NS)  clearly  separate  quarry  blasts  from  earthquakes.  A 
second  set  of  earthquakes  (the  cluster  located  within  the  network)  and  a  second  set  of  quarry  blasts 
from  the  same  mine  have  essentially  the  same  separation  (Figure  3). 

The  current  Automated  Time  Frequency  Discriminant  (now  known  as  the  ATFD)  recognizes  long  lived 
modulations  by  applying  three  separate  tests.  The  algorithm  estimates  the  time-independence  of  the 
binary  patterns  by  calculating  the  autocorrelation  of  individual  narrow  band  (single  frequency)  timeseries. 
The  “bandedness”  is  estimated  by  calculating  the  coda  cepstrum  and  taking  the  extreme  value  which  is 
independent  of  time.  Independence  from  recording  direction  is  judged  by  cross-correlation.  These 
individual  parameters  are  merged  into  a  single  discrimination  score  with  the  aid  of  multivariate  statis¬ 
tics  (Seber,  1984).  The  means  by  which  the  binary  sonograms  are  reduced  to  a  few  discrimination 
parameters  is  reviewed  in  Appendix  A  of  this  report. 

Global  test  of  the  ATFD.  We  have  used  a  wide  variety  of  datasets  to  test  the  utility  of  the  ATFD. 
Specifically,  we  have  tested  the  adaptability  of  the  algorithm  to  various  styles  of  seismic  deployment 
(me/,  single  3-Component  stations,  single-  and  3-Component  networks  and  arrays)  varying  seismic 
instrumentation  (me/,  short  period  and  broadband),  different  geologic  settings  and  mining  practice.  In 
the  past  year  we  have  analyzed  the  following  datasets  (Figure  1): 

1)  single  (chemical)  explosions  and  quarry  blasts  recorded  by  a  sparse  3-Component  network 
(NRDC)  deployed  in  1987  on  a  high  Q  crust  in  central  Kazakhstan.  The  events  occurred 
within  300  km  of  the  stations.  Event  identification  was  based  on  C.  Thurber's  interpretation 
of  SPOT  photos  (Thurber  et  al,  1989)  and  the  fact  that  the  Kazakh  platform  is  seismically 
quiet. 

2&3)  earthquakes  and  quarry  blasts  recorded  by  a  dense  (10  stations  in  an  aperture  of  200 
km)  telemetered  3-Component  broadband  network  (KNET)  deployed  in  Kyrgyzstan  between 
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Figure  2.  The  10  station  Kyrgyz  telemetered  broadband  network  (grey  triangles).  Two 
clusters  of  earthqaukes  are  contrasted  with  two  separate  sets  of  quarry  blasts  which 
occurred  in  a  mine  located  just  northeast  of  the  station  CHM  (at  the  crosshairs).  The 
remote  set  of  events  (dataset  "A")  are  aftershocks  to  the  magnitude  7.4  Suusamyr 
thrust  earthquake.  The  second  set  of  earthquakes  (dataset  "B")  occurred  within  the 
network,  roughly  50  km  south  of  CHM. 

the  Kazakh  platform  to  the  north  and  the  seismically  active  Tien  Shan  to  the  south.  We  have 
analyzed  two  clusters  of  earthquakes  and  contrasted  them  with  two  sets  of  quarry  blasts  which 
occurred  at  the  same  mine  (Figure  2).  The  limestone  mine  was  located  by  Rob  Mellors 
(UCSD)  at  43.028°  N  and  74.888°  E.  Events  were  determined  to  be  quarry  blasts  on  the  basis 
of  their  location  and  the  origin  time  -  the  limestone  shots  typically  occurred  in  the  late  after¬ 
noon.  In  addition  some  of  these  events  produced  obvious  acoustic  signals  which  arrived  at 
CHM  roughly  37  s  after  P.  The  network  made  triggered  recordings  at  100  sps  and  continuous 
recordings  at  20  sps. 
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figure  i.  hull  network  averaged  three-way  crosscorrelations  calculated  using  the  events  in  dataset  A  (left  tigures) 
and  B  (right).  We  find  that  the  network  solutions  exhibit  little  crossover  between  the  two  event  types.  Using  the 
earthquakes  in  dataset  A  we  have  calculated  a  mean  and  +/-  1,2  and  3  standard  deviations.  These  statistics  are 
overlain  on  the  network  solutions  from  dataset  B  showing  that  the  cross-correlations  of  the  two  sets  of  events  are 
essentially  the  same. 


4)  recordings  of  quarry  blasts  and  earthquakes  made  by  the  NORESS  small  aperture  array  in 
southern  Norway.  NORESS  is  comprised  of  25  stations  of  GS13  seismometers  in  a  3  km 
aperture  on  Precambrian  or  Paleozoic  hard-rock.  Four  of  the  stations  have  3-Components. 
The  signals  are  digitized  at  40  sps  The  events  in  this  dataset  occurred  within  7.5°  of  the  array. 

5)  recordings  of  quarry  blasts  and  earthquakes  made  by  GERESS  (a  NORESS  style  array  in 
Germany).  These  recordings  were  obtained  from  the  European  Ground  Truth  Database  and 
were  contributed  by  J.  Wuster.  The  earthquakes  largely  occurred  in  swarms.  The  quarry 
blasts  occurred  in  Western  Bohemia  (Wuster,  1993). 

6)  recordings  of  quarry  blasts  (in  the  Newmont  Gold  Mine)  and  earthquakes  (including  events 
in  the  Rock  Valley  earthquake  sequence  near  NTS)  made  by  a  single  3-Component  station 
(ELK  in  the  LNN  network)  located  near  Carlin,  NV  roughly  80  km  from  the  mine.  This  hard 
rock  mine  (which  detonates  ripple-fired  shots  every  2  days  each  using  15  to  100  tons  of 
explosives)  has  been  monitored  by  LLNL  scientists  (inch  Drs.  Peter  Goldstein,  Bill  Moran 
and  Steve  Jarpe).  Mining  records  have  been  used  to  constrain  the  blast  patterns,  which  typi¬ 
cally  involve  delays  from  50  to  100  ms.  These  data  were  provided  by  LLNL. 

The  ATFD  tested  on  a  dense  3-Component  network:  To  illustrate  the  technique  we  show  the  results  of 
the  analysis  of  the  second  and  third  datasets  (collected  by  KNET;  Figures  1  and  2).  The  upper  two  and 
lower  left  panels  of  Figures  4  and  5  we  show  the  raw  output  from  the  ATFD  (refer  to  appendix  A  for 
details).  These  panels  show  how  well  the  two  event  types  (earthquakes  and  quarry  blasts)  are  separated 
by  the  metrics  applied  to  the  time-frequency  expansions.  The  lower  right  panel  shows  the  linear  dis¬ 
criminant  scores  obtained  by  a  multivariate  analysis  of  the  raw  parameters.  Considering  dataset  A  we 
found  that  two  events  (of  47)  were  misclassified.  This  is  not  likely  due  to  poor  coverage,  since  one  of 
the  events  was  recorded  by  9  stations  (Figure  6),  but  likely  was  due  to  the  source  or  propagation.  As 
discussed  by  Hedlin  et  al.  ( 1989)  a  time  independence  can  be  acquired  through  propagation  by  reso¬ 
nance.  As  discussed  by  numerous  authors  (incl.  Stump  etal,  1994)  ripple  fired  shots  will,  sometimes, 
not  produce  scalloped  spectra.  Outliers,  such  as  these,  are  a  critical  problem  and  are  discussed  in  more 
detail  later  in  this  report. 

The  results  of  the  KNET  analysis,  and  the  other  studies  summarized  above,  are  shown  in  Figure  7 
where  we  show  the  misclassification  probabilities  (see  Appendix  B  for  a  full  description  of  these  analy¬ 
ses)  .  By  necessity  we  have  paid  the  closest  attention  to  the  problem  of  discriminating  quarry  blasts 
from  earthquakes.  We  expect  that  earthquakes  should  not  yield  coda  with  time-  or  component-inde¬ 
pendent  spectral  qualities  and  thus  have  used  them  to  serve  as  proxy s  for  single  explosions.  However, 
the  most  challenging  and  interesting  problem  involves  the  discrimination  of  single  and  multiple  explo¬ 
sions.  Althouth  we  haven’t  analyzed  a  large  number  of  single  explosions  (3  recorded  by  the  NRDC 
network  in  Kazakhstan;  dataset  #  1)  it  appears  that  the  ATFD  is  effective  at  separating  these  event  types 
(see  Hedlin  et  al.  ,1995  for  details).  To  further  address  this  problem  we  will  examine  the  single  explo¬ 
sions  recorded  during  the  1995  Black  Thunder  experiment  (Stump  et  al.,  1995;  Pearson  et  al,  1995). 

A  discriminant  based  on  a  wavelet  expansion:  In  the  preceeding  sections  we  analyzed  time  series 
using  multitaper  spectral  estimation  (Thomson,  1982).  using  data  selected  by  a  simple  sliding  window. 
While  this  approach  takes  advantage  of  multitaper's  proven  ability  to  yield  minimally  biased  spectral 
estimates  from  short  time  series  and  has  proven  useful  for  discrimination  it  has  several  drawbacks.  The 
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Figure  4.  KNET  A,  The  current  ATFD  yields  3  types  of  discrimination  parameters  -  all  of  which  are 
defined  in  Appendix  A.  One,  shown  in  the  lower  left  panel,  based  on  the  network  average  cross- 
correlation  across  the  three  recorded  channels,  is  similar  to  that  presented  in  the  preceding  figures.  A 
value  of  1  indicates  that  the  three  pairs  of  binary  sonograms  are  perfectly  correlated  at  all  stations. 
The  parameter  which  is  based  on  the  autocorrelation  of  single  frequency  timeseries,  is  shown  in  the 
upper  left  panel.  A  value  of  1.0  means  that  all  frequencies  at  all  stations  are  invariant  wrt  time  (they 
have  either +1  or  -1  at  all  times).  The  third,  based  on  the  coda  cepstrum  (see  Hedlin  et  ai,  1989;  1995), 
is  shown  in  the  upper  right.  A  binary  sonogram  which  is  periodic  in  frequency  and  independent  of  time 
would  have  an  extreme  value  of  1.0.  All  three  of  these  discriminants  show  good  separation  between 
the  two  types  of  events  in  dataset  A.  The  linear  discriminant  scores,  obtained  through  multivariate 
statistics,  are  shown  on  the  lower  right. 


resolution  in  time-frequency  expansions  based  on  sliding  window  ffts  is  invarient  with-respect-to  time 
and  frequency.  Typically  parameters  used  to  calculate  these  sonograms  are  chosen  in  an  ad  hoc  fash¬ 
ion. 

An  alternate  means  to  expand  a  time  series  into  a  time-frequency  plane  is  based  on  wavelets  {Daubechies, 
1990).  A  wavelet  transform  does  not  fix  the  tradeoff  between  time  and  frequency  localization  by  fixing 
the  time  window.  Instead  it  allows  the  time  and  frequency  resolution  to  vary  inversely  to  one  another 
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Figure  5.  The  same  as  Figure  4  except  this  figure  shows  results  from  the  events  in  KNET 
dataset  B. 


as  it  traverses  the  time  series.  The  wavelet  transform  has  been  shown  to  be  particularly  useful  in  the 
analysis  of  waveforms  beset  with  both  abrupt  changes  in  time  and  narrow,  highly-concentrated  fre¬ 
quency  bands  (similar  in  concept  to  the  tapers  used  in  multi-taper  analysis)  which  can  be  used  together 
to  yield  low-variance  time  and  frequency  representations  of  spectral  power  and  polarization.  Three 
wavelets  (calculated  assuming  a  time-bandcenter  product  of  3.0,  a  time-bandwidth  product  of  2.5,  an 
analysis  frequency  of  5.0  Hz)  are  shown  in  Figure  8. 

To  illustrate  the  potential  power  of  wavelets  we  consider  multitaper  and  multiwavelet  expansions  of  a 
controlled  signal.  We  use  the  recording  of  a  calibration  explosion  made  by  the  NRDC  station  BAY  (see 
Appendix  B)  as  an  empirical  Green’s  function  and,  assuming  linear  superposition,  synthesize  a  quarry 
blast  waveform  where  the  quarry  blast  consists  of  5  shots  spaced  at  50  ms  (Figure  9).  As  expected  the 
synthetic  sonogram  has  a  time-independent  spectral  high  at  20  Hz  (the  40  Hz  overtone  is  at  noise 
level).  In  Figure  10  we  show  the  multitaper  sonogram  and  time  series  of  16  seconds  including  the  P 
wave  onset  and  early  coda  from  the  calibration  explosion.  The  same  portion  of  the  synthetic  is  dis¬ 
played  in  the  middle  figure.  The  time-bandwidth  product  used  was  4  yielding  a  frequency  resolution  of 
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Figure  6.  Considering  KNET  dataset  A  events  we  plot  the  linear  discrimination  scores 
plotted  against  the  number  of  stations  used  (the  earthquake  scores  have  been  multiplied  by 
- 1  so  that  all  symbols  appearing  below  the  solid  line  -  at  zero  -  represent  misclassifications). 

We  find  that  one  earthquake  recorded  by  9  stations  failed  the  test.  A  quarry  blast  recorded 
by  3  stations  was  also  misclassified. 

4  Hz.  In  the  right  panel  (Figure  10)  we  display  a  multi  wavelet  transform  of  the  synthetic  quarry  blast 
time  series.  While  this  expansion  has  roughly  the  same  frequency  resolution  as  provided  by  the  multitaper 
approach  the  time  resolution  is  much  greater.  At  20  Hz  the  wavelet  spans  0.5  seconds. 

Outliers  and  Hybrid  monitoring :The  Black  Thunder  seismo/acoustic  experiment.  To  improve  the  sci¬ 
ence  of  small-event  discrimination  we  need  to  understand  better  why  some  events  are  misclassified  and 
we  need  to  explore  better  data  collection  and  processing  techniques.  To  investigate  the  physical  rela¬ 
tionship  between  the  ripple-fired  source  and  regionally  recorded  seismic  and  acoustic  wavefields  we 
have  formed  a  collaboration  with  Brian  Stump,  Craig  Pearson  and  Rod  Whitaker  at  LANL.  The  pri¬ 
mary  result  of  this  collaboration  is  an  experiment  which  will  occur  in  July  and  August  of  this  year.  As 
part  of  this  experiment  Brian  Stump  and  Craig  Pearson  from  LANL  will  deploy  video  equipment, 
accelerometers  and  acoustic  gauges  inside  the  Black  Thunder  coal  mine  in  the  Powder  River  Basin  of 
Wyoming.  The  Powder  River  Basin  includes  some  18  coal  mines  which,  in  total,  generate  an  average 
of  1  mine-blast  every  other  day.  The  in-mine  equipment  will  record  not  just  the  timing  of  the  shots  (t) 
in  the  mine  but  also  the  horizontal  spatial  positions  (x  and  y). 
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A  team  from  UCSD  will  deploy 
a  regional  network  of  broadband 
seismometers,  pressure  gauges 
and  infrasound  receivers  (Figure 
11).  The  network  will  include 
five  3-Component  broadband 
seismic  stations  (STS2’s)  sam¬ 
pling  continuously  at  100  sps. 
Four  of  the  stations  will  be  lo¬ 
cated  200  km  from  the  mine.  The 
fifth  station  will  be  deployed  100 
km  from  the  mine  on  the  same 
azimuth  as  one  of  the  outer  sta¬ 
tions.  This  temporary  network 
will  add  to  permanent  seismic 
deployments  that  include  the  3- 
Component  GSETT3  beta  station 
RSSD  and  the  array  at  Pinedale. 
Three  pressure  gauges  will  be 
deployed  at  the  stations  located 
200  km  to  the  north,  south  and 
west  of  the  mine.  Paired  with  the 
pressure  gauge  at  the  western  site 
will  be  a  small  infrasound  array 
consisting  of  3  infrasound  receiv¬ 
ers  in  a  triangle  100  m  on  a  side. 
In  part  we  will  be  pairing  the 
acoustic  and  infrasound  equipment  to  assess  the  signal  to  noise  gain  that  the  infrasound  receivers 
(loaned  to  us  by  Rod  Whitaker)  will  give  us.  The  acoustic  and  seismic  equipment  will  be  co-located  to 
allow  us  to  assess  the  value  of  a  seismo/acoustic  station.  The  acoustic  instrumentation  should  be 
located  within  the  first  acoustic  convergence  zone  and  we  expect  the  signals  from  the  mine  will  be 
relatively  strong  to  the  west  due  to  refraction  in  the  mid-atmosphere.  All  acoustic  signals  will  be 
digitized  at  40  sps. 

A  second  team  from  UCSD  will  deploy  3  continuous  GPS  receivers  to  detect  ionospheric  perturbations 
that  result  from  the  pressure  front  that  emanates  from  the  mine  and  is  not  ducted  between  the  mid¬ 
atmosphere  and  the  ground  but  propagates  into  the  ionosphere.  The  current  plan  has  2  receivers  lo¬ 
cated  near  to  each  other  at  a  range  of  50  km  to  the  southwest  of  the  mine.  A  third  receiver  will  be 
deployed  on  the  same  azimuth  but  at  a  range  of  200  km  (Figure  11). 

In  addition  to  this  near-regional  equipment  we  have  made  arrangements  with  Bob  Ellis  and  Michael 
Bostock  (both  at  the  University  of  British  Columbia)  to  collect  far-regional  recordings  of  the  mine 
blasts  at  700  km  (using  an  array  of  7  seismic  sensors  in  Saskatchewan,  Canada)  and  at  2000  km  (using 
an  array  of  seismic  sensors  located  in  the  Northwest  Territories,  Canada).  We  have  also  made  arrange¬ 
ments  with  Jim  Zollweg  (Idaho)  and  Scott  Smithson  (Wyoming)  for  additional,  near-  to  mid-regional 
recordings. 
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Figure  7.  Misclassification  probabilities  in  the  datasets  described 
in  Figure  1  and  reviewed  in  the  main  text  and  in  Appendix  B. 
Dataset  1  (NRDC)  contains  recordings  of  quarry  blasts  and  just 
3  single-explosions  and  is  not  summarized  statistically.  Datasets 
2  and  3  (KNET  A  and  B)  are  reviewed  in  detail  in  the  main  text. 


10 


best  3  wavelets  -  fc=5  Hz,  p=2.5,pc=3.0 


Figure  8.  The  three  best  bias  minimizing  wavelets  calculated  with  a  bandcenter  prod¬ 
uct  of  3.0;  a  time-bandwidth  product  of  2.5  and  an  analysis  frrequency  of  5.0  Hz. 


We  expect  that  this  instrumentation  will  record  four  significant  (cast)  explosions.  On  July  19  there  will 
be  a  shot  involving  ~4  million  pounds  of  explosives.  We  expect  two  additional  shots  (3  to  4  million 
pound  range)  to  occur  the  week  of  July  29  with  a  fourth  occuring  in  mid  August.  The  proximity  of  this 
mining  activity  with  regional  earthquakes  located  in  the  Intermountain  Seismic  Belt  to  the  west  of 
Pinedale  {Hansen,  1992)  should  provide  a  number  of  natural  event  records.  The  data  from  this  experi¬ 
ment  will  allow  us  to  test  the  utility  of  pressure  gauges  and  an  infrasound  array  to  detect  acoustic 
emissions  from  a  surface  mine  and  estimate  the  back  azimuth  to  the  source.  We  will  also  be  testing  the 
utility  of  the  hybrid  seismo/acoustic  monitoring  station  and  have  an  opportunity  to  address  the  issue  of 
outliers.  The  video  data  that  Brian  and  Craig  collect  will  enable  us  to  model  the  regional  waveforms 
and  should  allow  us  to  test  linear  superposition  at  four  azimuths.  Permanent  stations  and  the  temporary 
deployments  in  Canada  will  allow  us  to  determine  the  effective  range  of  the  ATFD.  In  total  we  hope  to 
be  able  to  answer  the  following  questions:  1)  Can  earthquakes  give  rise  to  time-independent  modula¬ 
tions?;  (2)  Under  what  circumstances  will  ripple-fired  events  not  yield  modulated  spectra?;  (3)  To  what 
extent  is  success  dependent  on  source-receiver  range  and  azimuth?  Using  linear  superposition  theory 
we  will  predict  the  time  independent  spectral  signature  that  should  be  present  in  mine  records.  Using 
our  multi-taper  sonogram  algorithm  {Hedlin  et  al,  1989)  we  will  expand  the  3-component  recordings 
into  time-frequency  displays  and  investigate  the  correspondence  between  the  expected  and  observed 
time  independent  signatures. 
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Chemex  2  recorded  at  BAY.CO 
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time  (s) 


Synthetic  quarry  blast  -  5x50  ms  delays 
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Figure  9.  A  sonogram  calculated  from  a  calibration  (single)  chemical  explosion  recorded  by  the  NRDC 
network  station  BAY  (see  Appendix  B)  is  shown  on  top.  On  the  bottom  we  show  the  sonogram  of  a 
synthetic  quarry  blast  which  comprises  5  shots  spaced  evenly  at  50  ms. 


sonogram  of  chemex  2  recorded  at  BAY  _ sonogram  of  synthetic  quarry  blast  -  5x.05  multiwavelet  time-varying  spectmm  -  p=2.pc=1 0.nwin= 
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W)[96  Regional  Deployment 


Figure  11.  The  1996  Black  Thunder  regional  deployment.  The  WY96  seismic  sta¬ 
tions  consist  of  3-Component  STS2  broadband  seismometers.  The  seismometers  will 
be  making  continuous  recordings  digitized  at  100  sps.  The  acoustic  sensors  are  DNA 
pressure  gauges.  The  W96  infrasound  array  will  be  deployed  200  km  to  the  west  of 
Black  Thunder  (at  MNTA).  This  array  will  consist  of  3  sensors  arranged  in  a  triangle. 

Software  contribution:  An  initial  version  of  our  ATFD  has  been  given  to  David  Jepsen  at  the  IDC  and  is 
currently  being  incorporated  into  their  Detection  Feature  Extraction  (DFX)  software  package. 

B.  CONCLUSIONS  AND  PLANS  FOR  THE  COMING  YEAR 


The  key  question  which  we  are  attempting  to  address  with  our  research  is:  Can  a  robust  discriminant  be 
developed  which  will  discriminate  large  numbers  of  industrial,  multiple,  explosions  from  nuclear  ex¬ 
plosions  and  earthquakes?  It  seems  that  the  generic  ATFD  is  robust  with  misclassification  probabilities 
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NORESS  B  -  a  test  of  a  vertical-  and  3-Component  array 


Figure  12.  The  NORESS  B  dataset.  This  datase,  provided  to  us  by  Jay  Pulli,  contains 
recordings  of  quarry  blasts  and  earthquakes  made  by  NORESS.  Some  events  are  located  as 
much  as  25°  from  the  array. 


ranging  from  0.5  to  3.5%.  The  ATFD  is  easily  adapted  to  a  wide  range  of  deployments  (from  single 
stations  to  arrays  and  networks)  and  thus  easily  automated.  Discrimination  using  time-frequency  ex¬ 
pansions  does  not  rely  on  expert  interpretation  but  is  routine,  robust  and  easily  automated. 

In  the  coming  year  we  will  broaden  the  scope  of  our  test  of  the  ATFD  by  analyzing  data  at  greater 
ranges  and  in  regions  not  closely  studied  before: 

•  Examine  data  collected  by  the  Saudi  Network  {Vernon  et  al,  1996) 

•  Examine  CNET  data 

•  Examine  J.  Pulli’s  NORESS  Ground  Truth  Database  to  test  ATFD  at  greater  range  (Figure  12). 

Small-event  discrimination  at  low  frequencies.  There  have  been  a  number  of  observations  of  a  spectral 
null  at  low  (<  10  Hz)  frequencies  {e.g.  Gitterman  and  Van  Eck,1993)  observing  mines  at  local  ranges  in 
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Israel).  A  null  is  predicted  by  linear  wavefield  superposition  and  is  due  to  source  finiteness.  If  this  is 
detectable  in  the  events  in  the  aforementioned  datasets  we  will  attempt  to  determine  if  this  spectral 
quality  can  be  used  for  discrimination  at  far-regional  distances  (O  400  km).  We  are  currently  analyzing 
NRDC  recordings  of  single  and  multiple  explosions  and  have  made  arrangements  to  analyze  record¬ 
ings  made  in  Wyoming  in  1995  by  LANL,  LLNL  and  AFTAC. 

We  will  examine  more  recordings  of  single  explosions  (e.g.  those  recorded  in  Wyoming  in  1995  by 
Stump  et  ai,  1995  and  Pearson  et  al,  1995). 

We  will  further  examine  wavelet  based  time-frequency  expansions  and  analyze  the  variations  of  ampli¬ 
tude  and  polarization  with  time  and  frequency  to  better  understand  the  genesis  of  seismic  coda.  We 
will  bring  a  wavelet  based  ATFD  on  line  if  it  proves  to  be  superior  to  the  existing  algorithm. 

We  will  continue  our  analysis  of  tightly  constrained  mining  events  and  dual  (seismic  and  acoustic) 
monitoring  technology  by  analyzing  the  data  collected  in  Wyoming  this  summer.  With  these  data  we 
will  address  the  question  of  the  physics  of  ripple-fired  wavefield  superposition,  the  range  to  which 
ripple-fired  events  can  be  detected  and  improve  our  understanding  of  outliers.  This  experiment  will 
enable  us  to  test  the  utility  of  continuous  GPS  for  event  detection. 
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Appendix  A  The  Automated  Time-Frequency  Discriminant: 

We  have  noted  (in  the  main  body  of  this  report  and  previously  in  Hedlin  etal,  1 989)  that  the  frequency 
content  of  seismic  onsets  and  coda  resulting  from  ripple-fired  explosions  is  often  highly  independent 
of  the  recording  component  and  time.  It  is  well  known  that  the  energy  is  often  scalloped  in  frequency 
(e.g.  Bell  &  Alexander,  1977;  Baumgardt  &  Ziegler,  1988).  These  qualities  are  due  to  source  finiteness, 
the  intershot  delays  or  a  combination  of  the  two  {Hedlin  et  al,  1989;  1990).  They  can  sometimes  be 
acquired  during  propagation  through  a  resonant  crust  {Hedlin  et  al.,  1989). 

To  determine  the  degree  to  which  these  qualities  are  present  in  a  seismic  coda  we  first  expand  the 
recording  into  a  time-frequency  display  (sonogram)  by  sliding  a  window  (usually  2.5  seconds  long) 
along  the  time  series  (with  the  window  sliding  20%  of  its  length  each  time).  A  spectral  estimate  is 
calculated  at  each  window  position  using  7  multitapers  {Thomson,  1982;  Park  et  al,  1987)  with  a  time- 
bandwidth  product  of  4.  Two  multitaper  sonograms  of  recordings  made  during  the  1987  NRDC  ex¬ 
periment  (Figure  Al)  are  displayed  in  Figures  A2  and  A3.  These  were  obtained  from  recordings  of  a 
single  (chemical)  explosion  and  a  quarry  blast.  As  described  in  Hedlin  et  al.  (1989)  each  multitaper 
spectral  estimate  is  converted  into  binary  form  through  convolution  with  two  boxcars  (usually  span¬ 
ning  4.4  and  2.0  Hz),  differencing  the  two  smoothed  spectra  and  replacing  all  locally  high  (and  low) 
spectral  values  with  +1  (and  -1).  Using  the  original,  spectral,  sonogram  we  estimate  the  average  pre¬ 
onset  noise  level,  N(f),  and  randomize  all  values  in  the  binary  sonogram,  B(f,t),  which  have  a  smaller 
spectral  amplitude  than  this  average.  This  randomization  is  recorded  by  a  mask,  R(f,t),  which  equals  1 
or  0  (indicating  a  randomized  point  or  one  that  is  untouched).  Noise  suppression  is  directed  at  long- 
lived  spectral  lines  which  would,  if  left  untouched,  give  a  sonogram  an  improperly  high  level  of  time- 
independence.  Four  binary  sonograms  are  displayed  in  Figure  A4.  The  upper  left  pattern  was  obtained 
from  a  recording  of  a  single  chemical  explosion  (Figure  A2).  The  other  three  patterns  were  obtained 
from  a  3-Component  recording  of  a  quarry  blast  (Figure  A3).  This  figure  shows  the  high  correlation  (in 
time  and  across  recording  directions)  of  the  scallops  present  in  the  binary  sonograms  extracted  from 
the  quarry  blast. 

Although  the  sonograms  contain  a  wealth  of  information  about  the  evolution  of  spectral  energy  in 
seismic  coda  they  must  be  collapsed  into  a  few  parameters  which  are  diagnostic  of  ripple-firing  and  can 
be  used  for  automatic  source  discrimination.: 

Independence  from  recording  direction.  This  quality  is  estimated  simply  by  calculating  the  zero  lag 
cross-correlation  between  the  three  pairs  of  binary  sonograms  (vertical  with  east-west;  vertical  with 
north-south,  north-south  with  east-west).  Given  two  binary  sonograms  E  and  Z  (obtained  from  the 
east-west  and  vertical  recordings  of  the  same  event)  the  zero  lag  cross-correlation  is  given  by; 

^  TTilr  E"ir  E{i,3)z{i,j)RE{i,mz{i,3) 

'ltt  re{%  j)Rz{ij) 

where  and  are  the  E-W  and  Z  randomization  masks.  This  sum  incorporates  only  those  points 
which  have  not  been  randomized.  The  perfect  correspondence  of  non-randomized  points  would  yield 
a  cross  correlation  of  1.0.  We  generally  do  not  rotate  the  data  into  physical  coordinates  since  this  was 
found  to  yield  no  improvement.  When  multiple  3 -Component  recordings  are  available  (either  in  a 
network  or  an  array)  the  individual  estimates  {e.g.  of  the  cross  correlation  between  Z  and  EW)  are 
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NRDC  -  A  sparse  3C  network 


Mines 

Calibration  explosions 


Figure  Al.  The  NRDC  network  was  deployed  in  1987  in  central  Kazakhstan  (see  Figure  1)  and  col¬ 
lected  triggered  250  sps  recordings  at  the  surface  and  in  boreholes  of  three  calibration  explosions, 
numerous  local  and  regional  quarry  blasts  and  a  large  number  of  far-regional/teleseismic  earthquakes. 
This  study  is  concerned  with  events  recorded  within  5°  of  the  stations. 

averaged  together.  Thus,  regardless  of  the  nature  of  the  seismic  3-Component  deployment,  component 
independence  yields  three  parameters  for  discrimination. 

Time  Independence  and  Periodicity  in  Frequency.  These  qualities  are  estimated  together  by  using  the 
coda  cepstrum.  The  cepstrum  is  calculated  by  taking  the  Fourier  transform  of  the  log  (base  10)  of  a 
single,  detrended,  spectral  estimate  {Baumgardt  &  Ziegler,  1988).  The  coda  cepstrum  (as  defined  by 
Hedlin  et  al,  1995)  is  calculated  by  taking  the  two-dimensional  Fourier  transform  of  the  binary  sonogram 
(which  is  detrended  by  the  conversion  into  binary  form).  The  two-dimensional  coda  cepstrum  is  sen¬ 
sitive  to  spectral  periodicities  which  are  constant,  or  cyclic,  with  time  in  the  onset  phases  and  coda.  To 
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Figure  A2.  Seismogram  resulting  from  a  single  chemical  explo¬ 
sion  (CH2;  Figure  Al)  detonated  in  Kazakhstan  and  correspond¬ 
ing  sonogram.  The  recording  was  made  at  a  range  of  157  km  by 
the  vertical  component  seismometer  at  Bayanaul. 


Figure  A3.  Seismogram  resulting  from  a  ripple-fired  quarry  blast 
(EVC;  Figure  Al)  detonated  in  Kazakhstan  and  corresponding 
sonogram.  The  recording  was  made  at  a  range  of  264  km  by  the 
vertical  component  seismometer  at  Bayanaul. 
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Figure  A4.  Binary  versions  of  the  sonograms  presented  in  Figure  B2  (upper  left)  and  B3  (upper  right).  Binary 
sonograms  calculated  from  the  north  lower  (left)  and  east  (lower  right)  component  recordings  of  event  c  (Fig¬ 
ures  B1  and  B3). 


estimate  the  degree  to  which  time-independent  scallops  are  present  we  take  the  maximum  of  the  coda 
cepstrum  at  a  time-frequency  of  zero.  This  value  is  normalized  relative  to  the  maximum  of  the  coda 
cepstrum  calculated  from  a  synthetic  binary  sonogram  (which  has  the  same  periodicity  in  frequency 
and  is  independent  of  time  with  the  exception  that  the  randomization  mask  applied  to  the  real  sonogram 
is  also  applied  to  the  synthetic).  The  quefrency  of  the  maximum  can  be  used  to  estimate  a  dominant 
ripple-fire  delay  or  source  duration.  When  the  event  has  been  recorded  by  numerous  sensors,  the 
cepstral  values  are  averaged  together.  Thus,  when  applied  to  a  single  component  deployment,  the  coda 
cepstrum  will  yield  a  single  parameter.  A  3-Component  deployment,  averaged  in  the  same  fashion, 
will  yield  three. 
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Figure  AS.  The  autocorrelations  calculated  for  three  quarry 
blasts  and  two  calibration  explosions  recorded  during  the  1987 
NRDC  experiment.  Also  plotted  are  the  autocorrelations  for 
synthetic  sonograms  that  are  perfectly  random  and  one  that 
is  independent  of  time. 


For  an  additional  measure  of  time-independence  we  apply  an  autocorrelation  operator  to  the  binary 
matrix,  where  the  autocorrelation  at  a  lag  of  k  windows,  A(k),  is  given  by: 

^  E-ir  B{i,j)B{i,j + k)R{i,j)R{i,j + k) 

Considering  events  recorded  by  the  NRDC  network  we  have  found  that  the  autocorrelation  does  not 
depend  on  the  lag  (provided  that  the  lags  are  large  enough  to  ensure  no  window  overlap).  In  practice, 
when  using  2.5  second  windows  with  20%  overlap,  we  use  minimum  lag  of  6  windows  and  a  maxi¬ 
mum  equal  to  the  duration  of  the  sonogram).  All  pairs  which  include  a  point  that  has  been  randomized, 
due  to  noise,  are  excluded.  For  a  robust  estimate  of  A  we  average  over  all  appropriate  lags.  As  above, 
autocorrelation  yields  a  single  averaged  parameter  (or  three  if  the  deployment  has  3-Component  sta¬ 
tions). 

With  the  measures  described  above  the  ATFD  yields  nine  parameters  when  applied  to  a  3-Component 
deployment  or  just  two  (from  the  coda-cepstrum  and  autocorrelation)  when  single  component  data  are 
available.  These  parameters  are  typically  merged  using  multivariate  statistics.  The  ATFD  is  trained  in 
a  new  region  by  processing  a  number  of  known  quarry  blasts  and  earthquakes  and/or  calibration  explo- 
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sions.  In  effect,  the  ATFD  is  taught  to  recognize  the  binary  patterns  produced  by  the  two  types  of 
events  (ripple-fired/non-ripple-fired). 

Appendix  B  —  A  Global  Test  of  the  ATFD 

We  briefly  described  the  results  of  the  KNET  experiment  in  the  main  body  of  this  report.  The  five  other 
studies  listed  in  section  two  are  briefly  summarized  in  this  appendix.  These  datasets  were  selected  so 
that  we  could  determine  if  the  ATFD  is  easily  adapted  to  different  styles  of  deployment  (single  3- 
Component  station;  tight  array;  sparse  or  dense  3-C  network)  and  successfully  discriminate  ripple- 
fired  events  from  single  explosions  and  earthquakes  in  spite  of  varying  mining  practice  and  geologic 
settings.  In  all  cases  standard  tests  were  applied  to  ensure  the  validity  of  the  statistical  analysis. 


24 


z{+).n(o).e(*)  z(+).n(o).en 


zn(+),ze(o),nen 


‘  0  2  4  6  8  10  12  14  16  18 

event  number 


Figure  Bl.  The  three  discrimination  parameters  (described  in  Appendix  A)  calculated 
using  the  NRDC  recordings.  The  quarry  blasts  are  events  1  through  15,  the  calibration 
explosions  are  events  16  to  18.  The  events,  with  the  exception  of  #s  1  and  10,  separate 
quite  well.  The  quarry  blasts  are  presumed  to  be  ripple-fired. 


NRDC:  The  1987  NRDC  experiment  was  unique  for  several  reasons.  Recordings  of  local  mining 
activity  and  teleseismic  natural  seismicity  were  made  in  Kazakhstan,  USSR,  by  a  sparse  3-Component 
network  situated  near  the  Semipalatinsk  nuclear  test  site  (Figure  Al).  Three  calibration  explosions 
(ranging  in  size  from  10  tons  to  20  tons)  were  detonated  during  the  test  at  regional  distances  from  the 
NRDC  stations.  A  disadvantage,  from  the  perspective  of  demonstrating  robustness  of  a  discriminant, 
is  that  just  3  single-shots  were  recorded.  We  can  demonstrate  visually  a  tendancy  for  the  ATFD  to  be 
able  to  separate  single  and  multiple  explosions  but  can't  support  this  with  a  statistical  analysis.  As  a 
result  we  display,  in  Figure  B 1,  the  raw  output  from  the  ATFD  but  no  multivariate  statistics.  What  we 
see  in  this  figure  is  a  clear  separation  of  most  of  the  quarry  blasts  from  the  calibration  events.  Two 
events  (1  and  10),  which  we  have  presumed  are  ripple-fired  quarry  blasts,  are  judged  to  be  similar  to 
the  calibration  events. 
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NORESS  A  -  a  test  of  a  vertical-  and  3-Comonent  array 


Figure  B2  Dataset  number  4.  Collected  by  the  NORESS  small  aperture  array.  This  dataset 
contains  recordings  of  26  quarry  blasts  (from  3  mines  including  Blasjo  and  Titania)  and  16 
earthquakes.  The  events  range  from  2.7  to  7.5°  from  the  array  and  have  magnitudes  of  1.6  to 
3.0.  A  preliminary  analysis  of  these  events  is  presented  in  Hedlin  et  al.  (1990). 


NORESS  A:  The  regional  NORESS  dataset  (first  analyzed  by  Hedlin  et  al.,  1990)  consists  of  42  events, 
26  quarry  blasts  and  16  earthquakes  (Figure  B2).  As  displayed  in  Figure  B3  all  statistical  measures 
applied  by  the  ATFD  show  an  obvious  separation  of  the  two  event  types  -  no  events  were  misclassified 
(the  misclassification  probability  is  3.7%).  This  is  a  promising  result  given  that  one  of  the  quarries  is 
located  7.5°  from  NORESS. 
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Figure  B3.  The  four  panels  in  this  figure  show  the  output  of  the  ATFD  applied  to  NORESS  recordings 
of  the  events  displayed  in  Figure  B6.  Events  1  through  26  are  quarry  blasts,  27  through  42  are  earth¬ 
quakes. 


VOGTLAND  -  a  test  of  a  vertical-  &  3-Component  array 


Figure  B4.  The  Vogtland  dataset  was  recorded  by  GERESS  and  obtained  through 
the  European  Ground  Truth  Database.  The  data  was  contributed  to  this  database  by 
Jan  Wuster. 


Vogtland:  The  events  in  this  dataset  were  recorded  by  GERESS  (Figure  B4).  They  were  obtained 
from  Jan  Wuster  via  the  European  Ground  Truth  Database  {Grant  &  Carabajal,  1995).  Included  are  1 1 
quarry  blasts  and  10  earthquakes.  All  events  are  located  between  1.5  and  2.5°  to  the  northwest  of  the 
array.  These  populations  are  easily  separated  by  the  ATFD  (misclassification  probability  of  0.55%; 
Figure  B5). 
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Figure  B6.  The  Carlin  dataset  includes  recordings  of  mine  blasts  made  at  the  Newmont  Gold  Mine  in 
Nevada  and  widely  distributed  earthquakes.  The  recordings  were  made  by  a  single  3-Component 
station  (ELK).  The  data  were  provided  to  us  by  Steve  Jarpe  and  Peter  Goldstein  at  LLNL. 

Carlin:  This  dataset,  contributed  to  us  by  Steve  Jarpe  and  Peter  Goldstein  at  LLNL,  consists  of  42 
quarry  blasts  (all  occuring  at  the  Newmont  Gold  Mine  in  Nevada)  and  36  widely  dispersed  earthquakes 
(Figure  B6).  The  events  occurred  within  2°  of  the  single  3-Component  station  ELK.  Again,  these 
events  are  easily  separated  by  the  ATFD.  Of  the  78  events,  just  1  was  misclassified  (the  misclassification 
probability  is  2.6%;  Figure  B7). 
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cariin  -  Linear  discriminant  function 
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PO  BOX  808,  MS  L-200 
LIVERMORE,  CA  94551 


LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN;  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  808,  MS  L-221 
LIVERMORE,  CA  9455 1 


DONALD  HELMBERGER 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 

DIVISION  OF  GEOLOGICAL  &  PLANETARY  SCIENCES 

SEISMOLOGICAL  LABORATORY 

PASADENA,  CA  9 1125 

ROBERT  HERRMANN 
ST.  LOUIS  UNIVERSITY 

DEPARTMENT  OF  EARTH  &  ATMOSPHERIC  SCIENCES 
3507  LACLEDE  AVENUE 
ST.  LOUIS,  MO  63103 

ANTHONY  lANNACCHIONE 
BUREAU  OF  MINES 
COCHRANE  MILL  ROAD 
POBOX  18070 

PITTSBURGH,  PA  15236-9986 
THOMAS  JORDAN 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
EARTH,  ATMOSPHERIC  &  PLANETARY  SCIENCES 
77  MASSACHUSETTS  AVENUE,  54-918 
CAMBRIDGE,  MA  02139 

LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  808,  MS  L-207 
LIVERMORE,  CA  94551 


LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

LLNL 

PO  BOX  808,  MS  L-175 
LIVERMORE,  CA  94551 
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LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  808,  MS  L-208 
LIVERMORE,  CA  94551 


LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  808,  MS  L-202 
LIVERMORE,  CA  94551 


LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  808.  MS  L-195 
LIVERMORE,  CA  94551 


LAWRENCE  LIVERMORE  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  808,  MS  L-205 
LIVERMORE,  CA  94551 


THORNE  LAY 

UNIVERSITY  OF  CALIFORNIA,  SANTA  CRUZ 
EARTH  SCIENCES  DEPARTMENT 
EARTH  &  MARINE  SCIENCE  BUILDING 
SANTA  CRUZ,  CA  95064 

DONALD  A.  LINGER 
DNA 

6801  TELEGRAPH  ROAD 
ALEXANDRIA,  VA  223 10 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  1663,  MS  F665 
LOS  ALAMOS,  NM  87545 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  1663,  MS  C335 
LOS  ALAMOS,  NM  87545 


KEITH  MCLAUGHLIN 
MAXWELL  TECHNOLOGIES 
P.O.  BOX  23558 
SAN  DIEGO,  CA  92123 


RICHARD  MORROW 

USACDA/IVI 

320  2 1ST  STREET,  N.W. 

WASHINGTON,  DC  20451 


JAMES  NI 

NEW  MEXICO  STATE  UNIVERSITY 
DEPARTMENT  OF  PHYSICS 
LAS  CRUCES,  NM  88003 


JOHN  ORCUTT 

INSTITUTE  OF  GEOPHYSICS  AND  PLANETARY  PHYSICS 
UNIVERSITY  OF  CALIFORNIA,  SAN  DIEGO 
LA  JOLLA,  CA  92093 


ANATOLI  L.  LEVSHIN 
DEPARTMENT  OF  PHYSICS 
UNIVERSITY  OF  COLORADO 
CAMPUS  BOX  390 
BOULDER.  CO  80309-0309 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  1663,  MS  F659 
LOS  ALAMOS,  NM  87545 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  1663,  MS  D460 
LOS  ALAMOS,  NM  87545 


GARY  MCCARTOR 

SOUTHERN  METHODIST  UNIVERSITY 
DEPARTMENT  OF  PHYSICS 
DALLAS.  TX  75275-0395 


BRIAN  MITCHELL 

DEPARTMENT  OF  EARTH  &  ATMOSPHERIC  SCIENCES 
ST.  LOUIS  UNIVERSITY 
3507  LACLEDE  AVENUE 
ST.  LOUIS,  MO  63 103 

JOHN  MURPHY 
MAXWELL  TECHNOLOGIES 
11800  SUNRISE  VALLEY  DRIVE  SUITE  1212 
RESTON,  VA  22091 


CHARLES  ODDENINO 
BUREAU  OF  MINES 
810  7THST.NW 
WASHINGTON,  DC  20241 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K6-48 
RICHLAND,  WA  99352 
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PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K7-34 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K6-40 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K7-22 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K5-72 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K6-84 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K5-I2 
RICHLAND,  WA  99352 


FRANK  PILOTTE 
HQ/AFTAC/TT 
1030  S.  HIGHWAY  AlA 
PATRICK  AFB,  FL  32925-3002 


JAY  PULLI 

RADIX  SYSTEMS,  INC. 
6  TAFT  COURT 
ROCKVILLE,  MD  20850 


DAVID  RUSSELL 
HQ  AFTAC/TTR 
1030  SOUTH  HIGHWAY  AlA 
PATRICK  AFB,  FL  32925-3002 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5704 

MS  0979,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0979 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5791 

MS  0567,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0567 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5704 

MS  0655,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0655 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  61 16 

MS  0750,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0750 


KEITH  PRIESTLEY 
DEPARTMENT  OF  EARTH  SCIENCES 
UNIVERSITY  OF  CAMBRIDGE 
MADINGLEY  RISE,  MADINGLEY  ROAD 
CAMBRIDGE,  CB3  OEZ  UK 

PAUL  RICHARDS 
COLUMBIA  UNIVERSITY 
LAMONT-DOHERTY  EARTH  OBSERVATORY 
PALISADES,  NY  10964 


CHANDAN  SAIKIA 

WOOODWARD-CLYDE  FEDERAL  SERVICES 
566  EL  DORADO  ST.,  SUITE  100 
PASADENA,  CA  91 101-2560 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  61 16 

MS  0750,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0750 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  93 11 

MS  1 159,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-1 159 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5736 

MS  0655,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0655 

THOMAS  SERENO  JR. 

SCIENCE  APPLICATIONS  INTERNATIONAL 
CORPORATION 
1 0260  CAMPUS  POINT  DRIVE 
SAN  DIEGO,  CA  92121 
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AVI  SHAPIRA 
SEISMOLOGY  DIVISION 

THE  INSTITUTE  FOR  PETROLEUM  RESEARCH  AND 
GEOPHYSICS 

P.O.B.  2286,  NOLON  58122  ISRAEL 

MATTHEW  SIBOL 
ENSCO,  INC. 

445  PINEDA  COURT 
MELBOURNE,  FL  32940 


JEFFRY  STEVENS 
MAXWELL  TECHNOLOGIES 
P.O.  BOX  23558 
SAN  DIEGO,  CA  92123 


DAVID  THOMAS 
ISEE 

29100  AURORA  ROAD 
CLEVELAND,  OH  44139 


LAWRENCE  TURNBULL 

ACIS 

DCI/ACIS 

WASHINGTON,  DC  20505 


FRANK  VERNON 

UNIVERSITY  OF  CALIFORNIA,  SAN  DIEGO 
SCRJPPS  INSTITUTION  OF  OCEANOGRAPHY  IGPP,  0225 
9500  GILMAN  DRIVE 
LA  JOLLA,  CA  92093-0225 

DANIEL  WEILL 
NSF 

EAR-785 

4201  WILSON  BLVD.,  ROOM  785 
ARLINGTON,  VA  22230 

RU  SHAN  WU 

UNIVERSITY  OF  CALIFORNIA  SANTA  CRUZ 
EARTH  SCIENCES  DEPT. 

1 156  HIGH  STREET 
SANTA  CRUZ,  CA  95064 

JAMES  E.  ZOLLWEG 
BOISE  STATE  UNIVERSITY 
GEOSCIENCES  DEPT. 

1910  UNIVERSITY  DRIVE 
BOISE,  ID  83725 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
8725  JOHN  J.  KINGMAN  ROAD 
FT  BEL  VOIR,  VA  22060-6218  (2  COPIES) 


ROBERT  SHUMWAY 
410  MRAK  HALL 
DIVISION  OF  STATISTICS 
UNIVERSITY  OF  CALIFORNIA 
DAVIS,  CA  95616-8671 

DAVID  SIMPSON 
IRIS 

1 6 1 6  N.  FORT  MEYER  DRIVE 
SUITE  1050 

ARLINGTON,  VA  22209 

BRIAN  SULLIVAN 
BOSTON  COLLEGE 
INSITUTE  FOR  SPACE  RESEARCH 
140  COMMONWEALTH  AVENUE 
CHESTNUT  HILL,  MA  02167 

NAFI TOKSOZ 

EARTH  RESOURCES  LABORATORY,  M.I.T. 
42  CARLTON  STREET,  E34-440 
CAMBRIDGE,  MA  02142 


GREG  VAN  DER  VINK 
IRIS 

1 6 1 6  N.  FORT  MEYER  DRIVE 
SUITE  1050 

ARLINGTON,  VA  22209 

TERRY  WALLACE 
UNIVERSITY  OF  ARIZONA 
DEPARTMENT  OF  GEOSCIENCES 
BUILDING  #77 
TUCSON,  AZ  85721 

JAMES  WHITCOMB 
NSF 

NSF/ISC  OPERATIONS/EAR-785 
4201  WILSON  BLVD.,  ROOM785 
ARLINGTON,  VA  22230 

JIAKANG  XIE 
COLUMBIA  UNIVERSITY 
LAMONT  DOHERTY  EARTH  OBSERVATORY 
ROUTE  9W 

PALISADES,  NY  10964 

OFFICE  OF  THE  SECRETARY  OF  DEFENSE 
DDR&E 

WASHINGTON,  DC  20330 


TACTEC 

BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 

COLUMBUS,  OH  43201  (FINAL  REPORT) 
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PHILLIPS  LABORATORY 
ATTN:  XPG 
29  RANDOLPH  ROAD 
HANSCOM  AFB,  MA  0173 1-3010 


PHILLIPS  LABORATORY 
ATTN:  TSML 
5  WRIGHT  STREET 
HANSCOM  AFB,  MA  01731-3004 


PHILLIPS  LABORATORY 
ATTN:  GPE 
29  RANDOLPH  ROAD 
HANSCOM  AFB,  MA  01731-3010 


PHILLIPS  LABORATORY 
ATTN:  PL/SUL 
3550  ABERDEEN  AVE  SE 
KIRTLAND,  NM  87117-5776  (2  COPIES) 
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